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We construct the chiral perturbation theory for two-color QCD with two quark flavors 
as an effective theory on the 50(6) / SO(5) coset space. This formulation turns out to be 
particularly useful for extracting the physical content of the theory when finite baryon 
and isospin chemical potentials are introduced, and Bose-Einstein condensation sets on. 
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1. Introduction 

The phase diagram of quantum chromodynamics (QCD) has attracted much at- 
tention in recent years. The region of high baryon density and low temperature is 
relevant for the description of deconfined quark matter, which can be found in the 
centers of neutron stars. It is expected to exhibit a variety of color-superconducting 
phases P 

Unfortunately, there is only very little firm knowledge concerning the behavior 
of the cold and dense quark matter. At very high densities, asymptotic freedom 
of QCD allows one to use weak-coupling methods to determine the structure of 
the ground state. On the other hand, the phenomenologically interesting region 
of densities corresponds to the strong-coupling regime where ab initio calculations 
within QCD are not available. 

At the same time, current techniques of lattice numerical computations are not 
able to reach sufficiently high densities, due to the complexity of the fermionic Dirac 
operator that occurs in the Euclidean path-integral measure. This gave rise to the 
interest in QCD-like theories that are amenable to lattice simulations, in particular 
the tw o-co lor QCD with fundamental quarks and three-color QCD with adjoint 
quarksPEl 

It turns out that these theories may also be studied by the means of a low-energy 
effective field theory similar to the chiral perturbation theory of QCD, and non- 
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trivial information about their phase diagram thus obtained. Within this approach, 
the structure of the phase diagram has been investigated in detail, including the 
effects of finite temperature! ^ ! '^ The model-independent predictions of the low- 
energy effective field theory have bee n complem ented by lattice computation 
and calculations within several models! '" '' I ' ' 

The aim of the present letter is to provide an alternative low-energy effective 
formulation of the simplest of this class of theories — the two-color QCD with two 
quark flavors. 51 While the general description of the whole class is based on the 
extended SU{2Nf) chiral symmetry of the underlying Lagrangian with Nf quark 
flavors, we construct the effective Lagrangian by exploiting the Lie algebra isomor- 
phism SU(4) ~ SO(6). We show that such a picture displays more transparently 
the physical content of the theory and at the same time allows for an easy deter- 
mination of the true ground state, which has been sought by a convenient ansatz 
previously. 

The paper is organized as follows. In the next section we summarize the basic 
features of two-color QCD to set the stage for the following considerations. Next 
we work out the mapping between the coset space SO(6)/SO(5) that we use, and 
the SU(4)/ 'Sp(i) used in the literature. The rest of the paper is devoted to the 
construction of the effective Lagrangian and its detailed analysis. 

2. Two-color QCD 

In this section we recall the basic properties of two-color QCD, following closely the 
treatment of Kogut et a/P^The distinguishing feature of two-color QCD is the pseu- 
doreality of the gauge group generators, the Pauli matrices, Tj* = —T^^Ti- Con- 
sider now a set of Nf quark flavors in the fundamental representation of the gauge 
group. As an immediate consequence, we may trade the right-handed component 
of the quark field, ipu (flavor and color indices are suppressed), for the left-handed 
conjugate spinor ip^ = oiTii\}* R . (The Pauli matrices <7fc act in the Dirac space.) 

Instead of the usual Dirac spinor, ip = (ip^ ipn) T , we now work with the left- 
handed spinor, ^> — (V>l ^'r) T , in terms of which the quark Euclidean Lagrangian 
of the massive two-color QCD at finite chemical potential becomes 



Here D v is the gauge-covariant derivative that includes the SU(2) gluon field. 
fl u is th e sta tic uniform external U(l) gauge field that incorporates the chemical 
potential. 1 ^ In the two-flavor case we shall deal with both the baryon number and 
the isospin chemical potential, /ib and /x/, respectively, so that f2„ will eventually 
be set to £l u = S v o(iibB + nil). Here, 



a The determinant of the Dirac operator is always real in two-color QCD, it is, however, positive 
only for an even number of flavors. 



C = i^a v (D v - fi„)W - to [i* T a 2 T 2 Af* + H.c] . 



(1) 





I 



l 
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are the baryon number and isospin generators, respectively. (The Pauli matrices Tfe 
act in the flavor space.) Finally 

-(ii) 

denotes the mass matrix in the basis of the spinor \l/ and cr„ stands for the four- 
vector of spin matrices, er„ = (— i, <Jk)- 

In the chiral limit and the absence of the chemical potential, the Lagrangian 
Eq. is invariant under the extended global symmetry SU(2Nf), which includes 
the naive chiral group SU(Nf)i, x SU(Nf)n and additional symmetry transforma- 
tions due to the pseudoreality of the gauge group generators. The global symmetry 
is spontaneously broken by the standard chiral condensate down to its Sp(2Nf) 
subgroup. 

The low-energy effective field theory for the Goldstone bosons of the broken 
symmetry is thus naturally constructed on the coset space SU (2N f) / Sp(2N ~f) . This 
is parametrized by an antisymmetric unimodular unitary matrix E, in terms of 
which the leading-order effective Lagrangian reads 

F 2 

£ off =^-Tr(V I ,EV„E t )-GReTr(JE). (3) 
The V's denote the covariant derivatives, 

v v E = ft / s-(n v E + sn2'), 

V !/ E t = d u Y$ + (Sto,, + OT£t), 

while J serves as a source field for E, and is eventually set to mM . The quark mass 
m is connected to the Goldstone boson mass squared m 2 by the Gell-Mann-Oakes- 
Renner relation 

mG = F 2 m\. 

It is worth emphasizing that the incorporation of the chemical potential into the 
effective theory involves no extra free parameters — the way the chemical potential 
enters the Lagrangian is fixed by the form of the covariant derivatives. 

3. The SO(6)/SO(5) coset space 

From now on we shall restrict our attention to the case Nf = 2. In that case, note 
the Lie algebra isomorphisms SU(4) ~ 5*0(6) and Sp(4) ~ 50(5). This allows us 
to recast the low-energy effective field theory on the 50(6)/50(5) coset P3I There 
are altogether five degrees of freedom, or Goldstone bosons, corresponding to the 
five independent entries of the antisymmetric unimodular unitary matrix E. 
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3.1. Matrix basis 

The mapping to the SO(6)/ 50(5) coset space is now provided by the formula 

E = niT,i, (4) 

where n is a six-dimensional real unit vector and Hj is a convenient set of indepen- 
dent antisymmetric 4x4 matrices. For E to be unitary, the basis matrices must 
satisfy the constraint 

EtE j + E]E i = 2^, (5) 
Such a relation is fulfilled for instance by the matrices 

*-':«)■ *-c?°)- Mr? 



E 4 



o ° ■ M"'?). M ° '? 
o -it 2 y V ir 2 o / V _lr 3 o 



This particular set has been chosen to comply with existing literature. In fact, 
Kogut et a/0 use the notation E c and E^ for our Ei and E 2 , respectively, while 
Splittorff et al^ denote our Ex, E 2 and E3 by Em, Eg and E/, respectively. 

Let us in addition show a simple argument that suggests how to choose in general 
a set of matrices satisfying Eq. Recall that six independent antisymmetric 
Hermitian 4x4 matrices generate the real Lie algebra 5*0(4) ~ 50(3) x 50(3). 
This means that we deal with two sets of three matrices, which can be shown to 
fulfill the usual anticommutator of Pauli matrices, {t^t^ } = 2Sij. By multiplying 
the matrices from one of the sets by i, we arrive at three Hermitian matrices, 
Hi = {E2, E3, Eg}, and three anti-Hcrmitian ones, Ai = {Ei, E4, E5}. These satisfy 
the relations 

{Hi ,Hj} = 2% , {Ai ,A j }= -2% , [Hi , Aj] = 0, 
that are equivalent to Eq. JSJ. 



3.2. Structure of the coset 

It remains to prove that Eq. J3J provides a one-to-one parametrization of the coset 
SU (4) / Sp(4) . To that end, note that any antisymmetric 4x4 matrix U may be 
expanded in the basis Ej, U — ZjEj, where Zi are in general complex coefficients. 
The unitarity of U constrains these coefficients as 

i Vftj 

the Xi and Xji being the real and imaginary parts of Zi, respectively. 

It is now crucial to observe that the products iEjEj for i ^ j span the set of 
15 linearly independent generators of 577(4) so that the unitarity of U requires 
separately \zi\ 2 — 1 and Xiyj — XjUi for all pairs of 



I 



I 
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£2, £4 


B = ±1;I = 


diquark and antidiquark 


£3, S5, e 6 


B = 0;I= ±1,0 


isospin triplet 7r 


Si 


B = 0;/ = 


singlet cr 



Fig. 1. Quantum numbers of the components of the matrix S in the expansion El , 

The latter condition means that the complex phases of all the Zi's must be equal 
so that Zi — riie lip with real rij, while the former one requires Yli n i = 1- ^ is 
a matter of simple algebra to calculate the determinant of U, 



detU = 




Since the elements of the coset SU(A)/Sp(A) are unimodular matrices, we are 
left with two distinct possibilities, ip = or <p = n/2. (The next solution, ip = tt, 
already corresponds to <p = with just the sign of all the n^s inverted.) 

In conclusion, every antisymmetric unimodular unitary matrix E may be cast 
in the form flf. where n is either real, or pure imagin ary yec t or . However, in the 
standard coset construction of the effective Lagrangianj " I " I the global symmetry 
group is required to act transitively on the parameter space of the Goldstone fields 
that is, the actual coset space must be connected. As the chiral condensate, above 
which we build our effective theory, is described by the matrix Si (note that M = 
—Si), we have to choose the connected component with real n, as in Eq. jJJ. 



3.3. Physical content of the basis matrices 

It is instructive to look at the transformation properties of the matrix E. This 
will allow us to classify the Goldstone modes by their baryon and isospin quantum 
numbers. 

Recall that E is an antisymmetric tensor under SU(4) that is, it transforms as 
E — » UY,U T for U G SU(A). For an infinitesimal transformation generated by the 
baryon number or the third component of the isospin we get 

(5 e E = ie(QE + EQ T ) = ie{Q, E}, Q = B,I. 

Let E be a general block matrix of the form ( ^ ^ ) . Then 

K \ [T ^ _ ( \{t^K} i[r 3 ,L] 



{B ^ } -{0 -N)> V-|[T3,M]-i{T3,iV-} 

The quantum numbers of the particular components of E are summarized in Fig. 
□ 

To gain more insight into the nature of the effective field E, let us assign to it 
a composite field, 

E -> i* T cT 2 T2E* + H.c, 
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which corresponds to the form of the mass term in Eq. Q). Such a composite 
operator may be regarded as an interpolating field for the Goldstone boson. 

With the explicit knowledge of the matrices Sj it is now straightforward to find 
the particle content of the corresponding interpolating fields, cf. also Fig. Q 

S 2 -> -^ T C l5 T 2 r 2 ^ + H.c, £ 4 -» -^ T C 15 T 2 t 2 ^ + H.c, 
4. Chiral perturbation theory 

We are now ready to write down the leading-order effective Lagrangian and use it 
to analyze the phase diagram of the theory. First, we have to minimize the static 
part of the Lagrangian in order to determine the ground state at nonzero chemical 
potential. 

4.1. Global minimum of the static Lagrangian 

From Eq. J3J) we can immediately infer the static part, 

Atat =-^-Tr[(0 [/ S + SO^)(Stfi 1/ +nJst)] -GReTr(JS). (6) 
We include the external source J in the general form 

J = jM, 

with real jj. Note that setting j\ — m, we reproduce the quark mass contribution 
to the effective Lagrangian. 

The other sources can be taken as infinitesimally small, since they essentially 
serve to generate the ground-state condensates, 

(Si> = dJT' 

From the orthogonality property, Tr(EjSj) = 4<5y, we find ReTr(JS) = 4j ■ n so 
that we have 

(S) = AGn. 

It is obvious that the vacuum condensate rotates on a sphere in the six- 
dimensional space, with coordinates corresponding to the six basis matrices Ej. 
It remains to calculate the vector n minimizing the static Lagrangian ©. 

Note first that in the absence of chemical potential, = 0, the static Lagrangian 
is minimal when the condensate is aligned with the external source j. 

When baryon and isospin chemical potentials are switched on, we shall for sim- 
plicity assume that only the sources j\ , j 2 , j'3 are present. This is sufficient to include 



I 



I 
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the quark mass effects and calculate both the diquark and the isospin (pion) con- 
densate. Taking into account the explicit form of the charge matrices, Eq. J3J, the 
static Lagrangian becomes 

£ stat = -2F 2 [n 2 B (n 2 2 + n\) + [ij(n 2 + n 2 )] - 4G(ji»i + Jan 2 + Jang). 

The first term is invariant under SO(2) x SO(2) rotations in the planes (2,4) 
and (3,5). In the absence of the external sources, this symmetry may be exploited 
to set 714 = "5 = 0. The source J breaks the symmetry and, in fact, prefers the 
solutions with 714 — n§ = 0. The problem of finding the ground state thus reduces 
to minimizing the expression, 

£ stat = -2F 2 (/j|«2 + A*/»!) - 4G(jwi + j 2 n 2 + j 3 n 3 ), (7) 

on the sphere S : n 2 = 1. 

The Lagrangian now does not depend on 714,715, so that we are actually 
looking for a minimum on the ball, ti 2 + n 2 . + 713 < 1. It is clear that at the global 
minimum, both terms on the right hand side of Eq. JJJ) are negative (otherwise we 
could lower the energy by the inversion, n — > —n). By the same token, the global 
minimum must lie on the surface of the ball, since if this were not the case, we could 
lower the energy by scaling up the vector: n — > tn, t > 1. 

We have thus shown that in the global minimum, n\ + n| + n| = 1 and 714 = 
n 5 = n 6 = 0, and the ground-state condensate is given by the linear combination 
E = riiSi + 112X2 + ^3X3. We stress the simplicity of the proof of this fact within 
the SO(6)/SO(5) coset formulation of the chiral perturbation theory. Indeed, using 
the standard SU (4) / Sp(A) formalism, Splittorff et al^ only assumed such a form 
of S, and also did not prove that the minimum thus found was global. 

To demonstrate the power of the formalism we have built so far, we shall next 
rederive the results of Kogut et aL^for the case of nonzero baryon chemical potential 
fiB- We shall thus set /ij = and ji = Sum. Isospin chemical potential can be 
introduced along the same lines and the results of Splittorff et al. would be easily 
recovered. 

With the assumptions made, the static Lagrangian becomes 

£ stat = -2F 2 ml(x 2 sin 2 a + 2 cos a), 

where x = /Jb/ttItt and a parametrizes the minimum, S = Si cos a + £2 sin a. (The 
same argument as above tells us that when \xj = and J3 = 0, then 713 = in the 
global minimum.) 

Now when x < 1, the minimum occurs at a — — only the chiral condensate is 
nonzero, this is the normal phase. When, on the other hand, x > 1, the Lagrangian 
is minimized by cos a = l/x 2 . In this case, the chiral condensate rotates into the 
diquark condensate as the angle a increases. This is the Bose-Einstein condensation 
phase. 
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4.2. Excitation spectrum 

The spectrum of excitations above the ground state is determined by the bilinear 
part of the Lagrangian. Expanding Eq. © in terms of the components rij , it acquires 
the form 

£ off = 2F 2 {d v n) 2 + AiF 2 p B (n 2 d n i - n 4 d n 2 ) - 2F 2 Li 2 B (n 2 + n\) - 4F 2 m 2 r m. (8) 
To proceed, we have to deal separately with the two phases of the theory. 



4.2.1. The normal phase 

When x < 1, the ground state expectation values of n are ri\ = 1 and all other 
components zero. The independent excitations above the ground state may be iden- 
tified with m, i = 2, . . . , 6, while n\ is expressed in terms of them via the constraint 
ft 2 = 1, 



"l 



\ 



1 — nf = 1 — 5 rt 2 + higher order terms. 



The bilinear part of the Lagrangian (JSJ) becomes 
^ = ^ (^n,) 2 + (cW - Li B N)(d N^ + „ fl tft) + V7V • ViVt + m 2 £ - 

i=3,5,6 



where we have introduced N = n 2 + iri4, a complex field that carries baryon number 
one. This field thus corresponds to the diquark, while N> describes the antidiquark. 
We find the following dispersion relations, 

E = \J p 2 + m\ pion triplet n%,n^,n§, 



E = ^/ p 2 + m 2 — [ib diquark N, 



E = yjp 1 + m 2 . + fiB antidiquark . 



4.2.2. The Bose-Einstein condensation phase 



For x > 1, the chiral condensate alone is no longer the proper ground state and the 
Bose-Einstein condensation sets. We therefore parametrize the field n as 

n = (pcos <p, p sin ip, 713, n 4 , n$, uq). 

The ground state corresponds to p = 1, (p = a and n 3 = n± = n 5 = n 6 = 0. 
We set (p = a + 9 so that the five independent degrees of freedom are now 9 and 
rij, i = 3, . . . , 6. The radial parameter p is given by 



\ 



1 — n 2 = 1 — i \_. n i + higher order terms. 



2=3 



I 



I 
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Fig. 2. The orientation of 8 is perpendicular to the direction of the ground-state condensate, 
which is represented by the dotted line. The coordinates ni,ri2 are labeled schematically by the 
chiral and the diquark condensate, respectively. 



The bilinear Lagrangian reads in this case, 



•Cbili 



2F 2 



= 2j(9 ! yni) 2 + (9 l /0) 2 +2i^B(^9on4-7i49o6') cosa+n B 

=3 



E ■ 

i— 3,5,6 



9 2 sin 2 a 



Three of the degrees of freedom, 773,775,776, again represent the pion triplet, now 
with the dispersion relation E — y/jp + jij^. The dispersions of the remaining two 
excitations are obtained by a diagonalization of the inverse propagator in the ((9, 774) 
sector. The result is 



Ms, 



El = f + + 3 cos^ a) ± ^ ^ fi 2 B (l + 3 cos 2 a) 2 + 160 2 cos 2 a, 



in accord with previous work 

+ 3 cos 2 a) 



The masses of these modes are given by 
3m 2 



2 
Ms 



= 0. 



In contrast to the normal phase, there is always one truly massless Goldstone 
boson stemming from the exact baryon number U(l) symmetry, which is sponta- 
neously broken by the diquark condensate £2. 

It is, however, worth emphasizing that the nature of this Goldstone boson, as 
well as of the massive mode, changes as the chemical potential increases. There are 
two reasons — the rotation of the ground state in the (ni, 772) plane, and the balance 
between the mass term in the bilinear Lagrangian and the term with a single time 
derivative. 

In the limit a — > 0, the parameter 9 is, to the lowest order, equal to 772 and the 
parametrization of the case x < 1 is recovered. Here the Goldstone boson is the 
diquark N = 8 + 7774. 

As the angle a grows, the orientation of 6 also changes in the (771,772) plane so 
that it is always perpendicular to the direction of the condensate, see Fig.0 In the 
limit a — > 7r/2 that is, [Lb ^> Tn w , the condensate is purely diquark and 8 has the 
quantum numbers of m, i.e. the a field. The Goldstone boson is now 774. Note also 
that it is a linear combination of the diquark and the antidiquark so that it has no 
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definite baryon number. This is, of course, not surprising since the baryon number 
is spontaneously broken and thus it cannot be used to label the physical states. 

5. Conclusions 

We have constructed the chiral perturbation theory for two-color QCD with two 
quark flavors on the SO(6)/SO(5) coset. We have provided an explicit mapping 
between this formulation and that used previously in literature, based on the 
SU(4)/Sp(4) coset space. 

The virtue of the present approach is that the orthogonal rotations, in contrast 
to the unitary symplectic transformations, can be easily visualized and the physical 
content of the theory thus made manifest. We were also able to give a simple proof 
of the fact that the condensate taken previously as an ansatz is indeed the true 
ground state, and we thus justified the assumptions made in the older work. 

Since the 5 l O(iV)-symmetric nonlinear sigma model is known to great detail, 
the connection provided in this letter can hopefully lead to the improvement in the 
understanding of the phase diagram of the two-color QCD at low energies, at least 
in the two-flavor case. 
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